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Abstract 

We seek to utilize the nonextensive statistics to the microscopic model- 
ing of the interacting many-investor dynamics that drive the price changes 
in a market. The statistics of price changes are known to be fit well by 
the Students-T and power-law distributions of the nonextensive statis- 
tics. We therefore derive models of interacting investors that are based 
on the nonextensive statistics and which describe the excess demand and 
formation of price. 

1 Introduction 

As is known from economics, the price of a security (for example) can 
be related to the law of supply and demand. That is to say, the excess 
demand is proportional to the price such that we can write approximately 

(1) 

and here A is the market depth [4l[5]. In the past decade or so, there 
have been many models [H [S] proposed that attempt to capture the dy- 
namics and statistics of market participants. These range from minority 
game models JT^, multi-agent models, and lattice super-spin models that 
encode the many degrees of freedom available to an interacting investor 
as the degrees of freedom of the variables and spins of the models. These 
models attempt to quantify the excess demand brought about by the mis- 
match of supply and demand between interacting investors in a market. 
The hallmark for the success of a model has been the ability of the model 
in reproducing the stylized facts of real markets. These are the heavy tails 
(power-law) of the distributions, anomalous (super) diffusion, and there- 
fore statistical dependence (long-range correlations) of subsequent price 
changes. 



Recently we reported on an application of the C.Tsallis nonexten- 
sive statistics to the S&P500 stock index [I] [3] . There we argued that 
the statistics are applicable to a broad range of markets and exchanges 
where anamolous (super) diffusion and 'heavy' tails of the distribution are 
present, as they are in the S&P500 ^Sj. In effect we have characterized 
the statistics of the price changes (the left hand side of Eq.(l) ) as being 
well-modeled by the non-extensive statistics. We now seek to examine 
the demand-side of the equation in light of our recent findings that the 
non-extensive statistics models well the statistics of the price changes in 
real markets. As such, we will seek to outline a method by which one can 
obtain many-investor models within the context of the Tsallis nonexten- 
sive statistics. We will derive our specific models utilizing the well-known 
techniques of the maximum entropy approach [B] [8] . Let us briefly review 
the maximum entropy approach to be utilized here. The nonextensive, 
least-biased probability distribution function (PDF) P{z, t) of an under- 
lying observable z{t) is obtained by maximizing an incomplete information 
theoretic measure equivalent to the Tsallis entropy [T] |6l E] 

(2) 

Here P{z, t) is the probability distribution function and will be shown to 
be of a power-law form, and is the degree of non-extensivity or equiva- 
lently the incompleteness of the information measure. The inverse of the 
normalization is the partition function Z{t), and f3(t) is a Lagrange mul- 
tiplier associated with the constraint(s). In order to build our model(s), 
we must specify the constraints to be utilized in our maximization pro- 
cedure. These constraints will be the known (or assumed) observables of 
interest, and that are presumed to capture the deterministic behavior of 
our many-investor system. An interesting model for the investors is the 
model of investor bias and demand developed by Cont and Bouchaud and 
generalized to super-spins by Chowdhury [HOIS]. This model assumes 
that the generalized spins are representative of the magnitude and direc- 
tion of demand (the bias) of the investors. We will adopt this model, 
with suitable modifications, as a first approximation for characterizing 
the deterministic investor dynamics as observables from within the con- 
text of the nonextensive statistics. Following the work of Chowdhury et. 
al. let us define the demand function of a system of interacting investors 
as a classical Hamiltonian-like function in physics. This Hamiltonian will 
then be dependent on the magnitude of the demand, and the other de- 
gree of freedom, the direction of the demand or the bias. The bias in 
this model will be taken in this initial model to be discrete and will be 
represented as a spin. Also, for simplicity, let us initially assume that the 
magnitude is fixed. The N-investor interaction potential in the Hamil- 
tonian then is taken to be made of discrete terms with simple constant 
ferromagnetic coupling strengths [Jij — J > 0) and an anti-ferromagnetic 
coupling {Lij — L > 0) which we treat in the global mean field sense of 
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the Bornholdt model (9] [10] such that 



(M = Nm., m, = ^) 

N N 
N N 

= — ,7 y ^ rii (Tj + i ^ ^ f i M 

(3) 

The first potential term can be seen to be the usual form of a ferro- 
magnetic spin-spin interaction with a coupling strength and models the 
herd-like, or collective opinion formation of investors. The spins represent 
discrete bias, though the generalization to a continuous degree of bias is 
straightforward, and can take on the values 

r +h Bull ( t) 

a, = i 0, Neutral (0) 
[ -1, Bear ( i ) 

and therefore when two investors minimize their risk and both agree in 
their bias aiaj = 1 the overall energy is minimized. Note that in this 
simple model if one of two investors is neutral, the interaction term is 
zero. The second potential term is anti-ferromagnetic (L > 0) and models 
contrarian investor behavior as in the Bornholdt model. We can simplify 
the interactions by assuming a mean field approximation. Now let us 
assume that the bias fluctuates. This is a reasonable assumption, and 
leads us to the necessity of describing the averages of the fluctuations in 
terms of the mean and the variance. These moments can be written for 
the ith investor as q-parametrized averages [B] [S| 

(5) 

In order to obtain a tractable solution that retains the essential interac- 
tion dynamics to lowest order, let us make a mean field approximation to 
our Hamiltonian. This will allow us to linearize the spin-spin interaction 
and will allow us to derive a closed form solution which, it is supposed, 
will allow us to qualify and approximately quantify the inherent behavior 
of the system. The full solution to this model would perforce involve the 
complicated spin-spin terms and will result in corrections to our solved 
form mean field solution. As stated then, the mean field solution will 
give us approximately the first order response of the system, given the 
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interactions. We then have {Hmf is the mean field Hamiltonian) 

jV jV 

i,j i 
N 

= -{J-L)M^au 
i 

m = and M = Nm 

N 



(6) 



Here m is the average bias per investor and can be seen to be the analog to 
the average magnetization per particle. Next let us build into this model 
the open nature of a market. That is, the number of investors is not a 
constant in a market, and we will account for this fact by the inclusion of 
the number of investors as a further constraint (observable) in our model. 
As such the observable (per investor) to be included in our maximization 
procedure will be 

K = - iL, 

and is a total investor number constraint multiplier which from 
the usual thermodynamic analogies goes as the 'chemical potential'. We 
therefore can write for the ith investor the following observables to be 
included in our entropy maximization (and dropping the i sub- and su- 
perscripts) 

(9) 

The maximum entropy approach then allows us to vary the entropy 
given the constraints such that 

^ ^ (10) 

and we obtain our least biased probability density function as a Tsallis 
non- extensive statistics power-law form 



p(^^ ^ {1 + ;%-!)[-(,/ -£)m. ,7 -il]}^ 



(11) 
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and here the partition function is related to the normalization and is given 
byZ, = ^aP(a) 

We now wish to examine the average bias in this model. Following 
the usual magnetic systems argument, the average bias can be written as 
(recall t= +1,4-= -1) 



m = - rij. = , 



(12) 



The question then is how to obtain (TV N 4) given that = Nj^ + 
No + A'^- . We can write as before 

rii = (rr) 

= P(a = +1) - P(iT = -1) , ^ 

(13) 

This expression can then be related to the average price change by the 
market depth and we obtain our desired result. That is, the Tsallis power- 
law statistical distribution for the price and price changes, as reported 
elsewhere [1] for stock market indices such as the S&P500 high frequency 
price data, is obtained from the individual investor bias distributions with 
a proportionality factor of market depth converting the excess demand to 
price in currency. We make use of the market depth A and write {x is the 
price) 

^ (14) 



2 continuous spin model 

We now wish to generalize beyond the limitations of the approximations 
we have built into the pure spin model. To do this, let us assume that 
the state vector for the system is of two dimensions, the magnitude and 
the bias. We will then work with continuum spins as in the Kosterlitz- 
Thouless Hamiltonian. Let us write down the observables for the inter- 
acting investors in the two dimensions (y= magnitude, &=bias 'angle') of 
magnitude and bias as V. We then have 



i .i i J 

(15) 

and here Vint is the Kosterlitz-Thouless interaction, relegated to the 
role of the potential. This total Hamiltonian with the inclusion of the 
previously discussed total number N of investors as a further observable 
will comprise the constraints in the maximization of the entropy that we 
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will perform next. But first, let us simplify the interaction term in the 
potential again by averaging over the jih spins such that (tt < 5 < here) 



A' 



= — ^ 2^ (^')- 

' (16) 

again we can maximize the non-extensive entropy (per-investor) given 
the constraints of the potential and the first two central moments of the 
fluctuating variables of individual magnitude and bias such that for the 
ith investor we have 



TV 

{{e-B))^^ = l'[i9-e)]P{9,tyd9 = o, 

[) 

TT 

{{0-0)')^ = j[{9-^'']P{B,tydB 



(17) 



The maximum entropy method then states that we must maximize 
the entropy given the observables as constraints. This yields the following 
variation of the g-averaged observables 



<5<5), - .5 {,S[(:..i„,„,>^ + {{B-B))^ + {{B - B)') ^ - |^(N)A ^ (1 



(18) 

The least biased PDF will again be of the non-extensive form and can 
be written as 

{1 + miQ - - + [- IT'"'" (^) + Ln^^om] - H)^ 

(19) 

where the partition function is again related to the normalization and is 
now Zq = J ndO 



We wish to obtain the average bias and relate it to the price change. 
We then write the regular statistical average as M and utilize the market 
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depth to obtain the correspondence 



M = {cosie))= co4e)Pid,t)de, 



TT 

j COS 



I) 



{Aa;(f)) 



(20) 



3 continuous model 



We now wish to relax all of the approximations and the simplifications we 
built into the previous spin based models. To do this let us assume that 
the state vector for the system is comprised of two continuous degrees of 
freedom, the magnitude and the bias. Without restating the problem, let 
us write down the Hamiltonian for the investors in the two dimensions 
(y=magnitude,6'=bias 'angle') of magnitude and bias as Ho (denotes the 
non-interacting investor Hamiltonian). Also, let us propose some general 
interaction potential V the form of which will be examined subsequently 
. The important point here is that we are seeking to cast the problem 
of building a model for many interacting-investors into the powerful lan- 
guage of the many-particle physics as we feel this allows us to map some 
questions of modeling financial markets and investor behavior directly to 
well known physics-based paradigms. We have already touched upon this 
in our specific models discussed above, and now we wish to generalized 
the application of the technique. 

The model will consist of the free and interacting parts of the Hamilto- 
nian. That is, the ith investor will have i{y{t), o-{t)) magnitude of demand 
, with a bias of buy, sell or hold, at time t. The magnitude of demand 
can be assumed to be the demand for number of shares of a stock, and 
with the total number being considered fixed as a long-term constant 
n- ~ y{t)dy. The direction of demand, is 

discrete (—1,0,-1-1) corresponding to buy, sell or hold. The Hamilto- 
nian is 



2 

] 

with the constraints of the moments 



(21) 



() (1 

{{y^i -y,,f) = E J [y^i - y.^Y ■i'{y^i,f^i,t)"d.y>,i 

" *' (22) 

the maximization of the entropy then obtains the least biased distribution 
and with the normalization, ^ J P{ai,t)dya. = 1 to unity. The extrem- 
ization yields the least biased probability distribution function, P(ai,t) 
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which upon taking the mean field approximation ^ < aj >= M gives 



1 + /3(f)(7 - l)[(y„, - y,^ Y - a(t).IMai + j(t)K!Ti] 

(23) 

The expectation value for the demand is the expression < a >— 
J2a,i I '^iPiy<Ti^^i^i)dy<Ti, which explicitly is 

jV max 

i o 

rriKX niftx 

= J Piy+i,ai = +lj,)dy+i- J P{y i,(Ti = -lj)dy^i 
u 

(24) 

in terms of the number of shares demanded to be bought or sold or 
held for the moment at time t, the expectation value of excess demand is 

l\+i{t) + N i{t) + Nu{t)' 

(25) 

This excess demand is proportional to the price and with the market depth 
A as the proportionality factor we obtain the price and change in price if 
the definition of the variables of excess demand is relative y — x — x' 



(26) 



This model is perhaps the most detailed of the three discussed . It is 
also the most robust and points the way to the inclusion of interaction 
terms that describe interactions observed factors in real world markets. 
These could include the independent investor, institutional investor clus- 
tering, floor trader and outside trader time lags, non-constant market 
depths, more complex herding behavior etc. These can all be included 
as observables multiplied by Lagrange multipliers as constraints in the 
maximization, and the derived least biased distributions though perhaps 
complicated can be solved numerically. 



4 conclusion 

The models presented here are of increasing complexity. However the 
three models, discrete and continuous can be generalized further and can 
easily be applied. The numbers of shares being bought, sold or held at 
any moment determine the instantaneous excess demand. This in turn 
is related to the instantaneous price or price change if relative variables 
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are used, by the market depth proportionahty factor, here assumed con- 
stant. The statistical distributions of price changes and price have been 
shown to be well fitted by the students-T distribution and more recently 
the Tsallis nonex- tensive statistics distribution . The information theo- 
retic approaches taken here assume the nonextensive entropy as a starting 
point, and obtains not surprisingly a power-law distribution for the num- 
bers of shares to bought, sold or held. These are then summed and the 
time-dependent price change distributions obtained. A question to be an- 
swered in subsequent work is , what form of nonlinear interaction causes 
to arise a power-law distribution of price changes from a Gibbs-Boltzmann 
form of extensive entropy or information measure. How do nonlinear in- 
teractions modify the Gaussian distributions obtained from the extensive 
entropy. Also, how does the numerical simulation of the model, by Monte 
Carlo or stochastic trajectory methods, compare to actual market data 
price changes. Previous numerical simulations of power-law distributions 
have been shown by us to fit very well the high frequency stochastic time 
series data of price changes of the S&P500 and we have also applied the 
power-law statistics to generalized Black-Scholes equations of prices for 
options and derivatives secondary markets for the underlying stocks and 
stocks indices etc. We expect the present models especially the power-law 
distributed model to fit accurately the market price data , and given the 
parameters extracted from market data, the model can be a theory that 
begins to describe the statistical uncertainties, and on the average the 
microscopic interactions that occur between investors in a financial mar- 
ket and which then the theoretical model can subsequently be utilized for 
minimization of risk and the design of investment strategies in economics 
and finance. The author wishes to state that most of this manuscript was 
written in 20 Nov 2001 with M.D. Johnson and John Evans at the Uni- 
versity of central Florida Physics Department, Orlando, Florida. Recent 
research of the author's has refocused on this area of research, and the 
publication follows. 
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